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DETERMINATION OF PROPERTIES OF AN OPTICAL DEVICE 

BACKGROUND OF THE INVENTION 

The present invention relates to the determination of properties of an optical 
device under test, e.g. the determination of the group delay of the optical 
device. The group delay is a fundamental property of optical devices, such as 
single mode optical fibers or optical components such as Bragg gratings which 
devices are used in the optical transmission of information. 

As for the background of optical transmission of data it has to be said that the 
premier feature of optical fiber is its extremely low loss. This has made it the 
dominant transmission medium for long link lengths. The loss characteristics 
of fiber determines were optical communication is practical. At 1 550 nanometer 
single mode optical fiber has an attenuation of 0.2 dB/km. This allows fiber 
optic signals to be propagated through very long length of fiber without 
regeneration. Telecommunication systems use the 1300 and 1550 nm windows 
for lowest loss in the fiber. Since a telecommunication system must cover a 
very large distance, the aforementioned attenuation of the single strength in the 
fiber is of high importance. Therefore, the loss characteristics of optical fiber 
often limit the distance that a signal can propagate in the fiber. 

However, this is not always the case. In single mode fibers, chromatic 
dispersion can limit the distance over which fiber optic signals can propagate. 
Chromatic dispersion describes the fact that the speed of signal propagation 
in the fiber depends on the wavelength of the light. The consequence is that as 
the signal propagates through a long length of fiber, the edges of the waveform 
of the signal start to become more rounded. Eventually, the adjacent bits start 
to overlap in time causing the digital waveform to have poor readability. The 
amount of signal rounding depends beyond other parameters on the amount 
of chromatic dispersion in the cable. The problem of this pulse spreading is a 
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problem in today's fiber communication because the increasing bit rate of the 
state of the art fiber communication systems, which bit rate reaches numbers 
of up to 40 Gbit/s per channel, brings the chromatic dispersion in the range of 
the bit resolution of such a 40 Gbit/s communication system. 

To compensate chromatic dispersion in a fiber chirped Bragg gratings were 
developed. In these gratings different wavelengths of the pulse need different 
amounts of time traveling through the grating to compensate the dispersion 
caused by the fiber. The concept of this compensation is schematically 
illustrated in Fig. 1 . Fig. 1 shows the original pulse 1 that is broadened to a 
spreaded pulse 2 while traveling through a fiber 3. With the help of a circulator 
4 a chirped Bragg grating 5 is introduced in the path of the pulse 2. As 
described before the resulting pulse 6 has the original shape of pulse 1 again. 

However, to use such a component as a chirped Bragg grating it is necessary 
to have exact knowledge of the chromatic dispersion caused by such a grating. 
The measurement of chromatic dispersion is accomplished by analyzing the 
group delay through the fiber as a function of wavelength. To measure the 
group delay a wavelength tunable optical source such as a tunable laser is 
used to generate coherent light at different wavelengths. In the prior art the 
wavelength of the tunable laser is then incremented step by step and for each 
wavelength step a group delay is determined. Finally, the group delay details 
are used to calculate the chromatic dispersion coefficient. The disadvantage 
of measuring the group delays step by step is the time, which is needed to 
perform all the wavelength steps and the respective measurements. In other 
words this means it takes quite a long time to get data, which are so precise 
that this precision is sufficient for the needs of the telecommunication industry. 

SUMMARY OF THE INVENTION 

Therefore, it is an object of the invention to provide an improved determination 
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of a property of an optical device under test, e.g. the group delay of the optical 
device. The object is solved by the independent claims. 

The main advantage of the present invention is that it is possible to measure 
the group delay and therefore the chromatic dispersion of an optical device, 
5 such as a fiber or a chirped Bragg grating, in a short amount of time while still 
keeping the precision of the data high enough for the needs of the 
telecommunication industry. This goal is reached by the invention by tuning the 
frequency of the coherent light beam of the laser from a maximum to a 
minimum of a given frequency range in a given time interval, e.g. making a 

10 wavelength sweep through the given wavelength range. Because of the 
inventive correction of effects caused by a non-linearity in the tuning of the 
laser this sweep can be done without detrimental broadening of the resulting 
spectral width. Method and apparatus of the invention avoid the 
aforementioned problems of the prior art and provide for exact data well 

1 5 keeping the measurement time low. 

In a preferred embodiment of the invention there is performed a filtering of the 
Fourier transformed first signal with a high pass filter. It is further preferred to 
use a Hanning window as a shape for the high pass filter. Using this filter 
makes sure that a good elimination of not usable data is possible while still 
20 having important parts of data at the edges of the spectrum within the filter. 
Moreover, it is preferred to use a half of a Hanning window as a shape for the 
high-pass filter. 

In a further preferred embodiment of the invention this high pass filter is 
adapted to the precision and the shape of the resulting spectrum of the signal 
25 corrected for a non-linearity of the laser by making an interferometric signal out 
of the corrected first phase signal, Fourier transforming the interferometric 
signal to get a spectral signal, determining a fraction of the maximum of the 
spectral signal, determining the abscissas of the intersections of the ordinate 
-3 
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of the fraction which the curve of the spectral signal, determining the mean 
frequency f mean as the average of the abscissas, band pass filtering the spectral 
signal with a band pass filter having its center at the mean frequency and 
having a width greater than the width of the frequency range. The inventive 
5 concept of adaptive filtering of the signal corrected for a non-linearity of the 
laser is doubling the inventive success since by eliminating the non-linearity the 
resulting peak of the Fourier transformed spectrum is very sharp so that it can 
be used a filter width which is very much smaller than known filter width in the 
prior art. As can be seen in the following detailed description of the invention 
10 and the respective drawings this inventive concept is resulting in very smooth 
curves showing the group delay. 

g In a preferred embodiment the inventive method comprises the following steps: 

m - splitting the first initial light beam into a first light beam and a second light 

p beam, 

jL 15 - coupling the first light beam into the optical device under test, 

HJ - letting the second light beam travel a different path as the first light beam, 

^ - superimposing the first and the second light beam to produce interference 
between the first light beam and the second light beam in a resulting first 
superimposed light beam, 

20 - detecting as a first signal the power of the first superimposed light beam as 
a function of time when tuning the frequency of the coherent light beam from 
a minimum to a maximum of a given frequency range in a given time interval, 

- splitting the second initial light beam in a third light beam and a fourth light 
beam, 

25 - superimposing the third light beam and the fourth light beam after each light 

-4 



Agilent Technologies, Inc. 
Our Ref.: 20010487 

beam has traveled a different path, to produce interference between the third 
and the fourth light beam in a resulting second superimposed light beam, 

- detecting as a second signal the power of the resulting second superimposed 
light beam as a function of time when tuning the frequency of the coherent light 
beam from a maximum to a minimum of a given frequency range in a given 
time interval, 

- using the detected second signal for deriving a non-linearity information about 
a non-linearity in a tuning gradient of the frequency when tuning the frequency 
of the coherent light beam from the maximum to the minimum of the given 
frequency range, and 

- using the non-linearity information for correcting effects on the first signal 
caused by the non-linearity to get a corrected first signal. 

It is further preferred to flatten the first and/or the second signal before 
transforming it by approximating it with polynoms of at least second order. 

Additionally it is preferred to tune the frequency with a tuning velocity of about 
5 nm/s. 

Additionally it is preferred to tune the frequency within a tuning interval of about 
c/1400nm to c/1600nm, c = velocity of light. 

Additionally it is preferred to detect the power at a rate of about 100kHz. 

The term "coherent" in this application means that the coherence length of the 
light beam is bigger than the difference of the path of the first and second and 
third and fourth light beams, respectively. 

Other preferred embodiments are shown by the dependent claims. 
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It is clear that the invention can be partly or entirely embodied or supported by 
one or more suitable software programs, which can be stored on or otherwise 
provided by any kind of data carrier, and which might be executed in or by any 
suitable data processing unit. 



BRIEF DESCRIPTION OF THE DRAWINGS 

Other objects and many of the attendant advantages of the present invention 
will be readily appreciated and become better understood by reference to the 
following detailed description when considering in connection with the 
accompanied drawings. The components in the drawings are not necessarily 
to scale, emphasis instead being placed upon clearly illustrating the principles 
of the present invention. Features that are substantially or functionally equal or 
similar will be referred to with the same reference sign(s). 

Fig. 1 shows a chirped Bragg grating in a light path, 

Fig. 2 shows an embodiment of the invention, 

Fig. 3 shows a flow chart of the steps of an embodiment of the inventive 
method, 

Fig. 4 shows an illustration of the step 100 of Fig. 3, 

Fig. 5 shows an illustration of the step 106 of Fig. 3, 

Fig. 6 shows an illustration of the step 1 16 of Fig. 3, 

Fig. 7 shows two possible shapes for filters used in step 1 1 6 of Fig. 3, 

Fig. 8 shows an illustration of the step 100 of Fig. 3, and 
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Fig. 9 shows the group delay of the grating 1 2 of the example of Fig. 1 . 

DETAILED DESCRIPTION OF THE INVENTION 

Referring now in greater detail to the drawings, Fig. 2 shows a schematic 
illustration of a preferred embodiment of an apparatus 10 for interferometric 
5 determination of the group delay of a chirped Bragg grating 12 as a device 
under test. However, with the present invention it is possible to determine 
optical properties, e.g. the group delay ,of other devices, e.g. fibers or even air, 
also. 

The apparatus 10 comprises a tunable laser 14 as a signal source for a 
'¥- 10 coherent laser beam 16, wherein the laser 14 is tunable between 1400 nm and 
g at most 1 600nm. The tunable laser 14 can be tuned continuously within a given 
H} frequency range in a given time interval. The laser beam 16 is coupled into a 
SI first beam splitter 18, which splits the coherent laser beam 16 in a first initial 
beam 20 and a second initial beam 22. The first initial beam 20 is coupled into 
=1 1 5 a second beam splitter 24. With the help of the second beam splitter 24 the first 
initial beam 20 is split into a first beam 26 and a second beam 28. The first 
% beam 26 is traveling a different and longer path than the second beam 28, 
O since the first beam 26 is coupled via a circulator 52 into the grating 12. 
Reflected by the grating 12 the first beam 26 is coupled via the circulator 52 in 
20 the first path 32 of the first beam 26 again. Then the first beam 26 and the 
second beam 28 are united with a third beam splitter 34 by superimposing the 
first beam 26 and the second beam 28 to a first resulting beam 36. The first 
resulting beam 36 is coupled into a first photodiode 38. The photodiode 38 
transmits its output as a first signal to an analog/digital-converter (ADC) 40 
25 (which can be National Instruments AD-MIO-16DE-10). The ADC 40 is 
connected to a transient recorder 42 that can be a FAST TR1202 of Fast 
Comtec GmbH, Germany. The laser 14 can be a HP81680A of the applicant. 
The transient recorder 42 is coupled to a computer 44 as an evaluation unit of 
-7 
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apparatus 10. 

Second beam splitter 24, first beam 26, second beam 28 and third beam 
splitter 34 built up a first Mach-Zehnder-interferometer. 

Exiting the first beam splitter 18 is a second initial beam 22. A fourth beam 
5 splitter 46 splits the second initial beam 22 into a third beam 48 traveling a third 
path and a fourth beam 50. The third beam 48 is traveling a different and 
longer path than the fourth beam 50, wherein the path length of the third path 
differs from the path length of the fourth path by more than 1 m, preferred by 
less than 20 m, more preferred by less than 9 m, even more preferred between 

10 7 and 9 m. This different path is symbolized by some loops 30. The third 48 
and the fourth 50 beam are coupled into a fifth beam splitter 54 which units the 
third beam 48 and the fourth beam 50 to a superimposed second resulting 
beam 56. The second resulting beam 56 is then coupled into a second 
photodiode 58. The second photodiode 58 transmits its output as a second 

1 5 signal to an analog/digital converter ADC 60, similar to the ADC 40. The ADC 
60 is also connected to the transient receiver 42. Fourth beam splitter 46, third 
beam 48, fourth beam 50 and fifth beam splitter 54 built up a second Mach- 
Zehnder interferometer. 

The system for determination of the group delay of the optical component 12 
20 according to the present invention is based on the use of aforementioned 
Mach-Zehnder interferometers, both of which are supplied with light of the 
same light source 14. It is necessary to make the determination by 
interferometric measurements since the phase characteristics of the 
component 12 can not be detected directly. 

25 The calculation of the properties based on the detected first and second signals 
is as follows: 
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The laser 14 emits light 16 with field strength E dependent on the time t 
according to the following formula: 

E w (0=E 0 -expO -2^ \v mom {f)dn (3.1) 

0 

5 with v (3.2) 
2k dt 

into the Mach-Zehnder interferometers, u mom denoting the instantaneous 
frequency of the laser 14 . 

In the embodiment of Fig. 2 the frequency of the laser 14 is changed in a linear 
1 0 way so that the instantaneous frequency is 

^ om (0 = v 0 +/-/ (3.3) 

To calculate as a second signal the power of the resulting signal 56 of the 
second Mach-Zehnder interferometer detected with the second detector 58, 
1 5 one puts for the detected voltage U det that the light of the longer path is delayed 
by x = AL n/c when arriving at the fifth beam splitter 54. n is the refractive index 
of the medium, Al_ is the difference of the path length and c is the velocity of 
light. 

The detected voltage is determined by the two superimposed fields that travel 
20 through the different paths of the interferometer: 
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(3.4) 



Only regarding of the changing part of this signal and putting r =fu(t) dt the 
result is: 

U~ et (t)=2-cos[2x-(r mom (t + r)-r mom (t))]. 0-5) 

Assuming x « t and approximating with differential quotient one reaches: 



The second signal of the second interferometer is therefore ideally, i.e. with a 
linear rise in frequency of the laser 14 over time a cosine, the frequency of 
which is determined by the length AL and the sweep velocity of the wavelength 
over time. 

Because the laser 14 will produce not an ideal sweep over the given frequency 
range and will produce non linearities in this sweep one can expect a cosine 
which has a frequency which is changing overtime. With this second signal it 
is possible according to the present invention to correct the first signal of the 
first interferometer for this non-linearities in therefore getting a resulting 
corrected first signal as if the laser 14 had performed a linear sweep. 




2-cos[2tzt- V mom {t) -t] (3.7) 



verage over i 
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In the first interferometer in which the grating 12 is incorporated the group delay 
of this grating 12 is determined. For this purpose as the first signal the detected 
power of the first photodiode 38 is calculated in an analog way as explained 
above. However, in the following the first signal is only determined in its 
5 dependency of the optical frequency u, where as the time dependency is 
determined later on. 

For the detected signal at the photodiode 38 the following formula applies: 

U det (v) ~ | exp (-7 • j3 DUT (v) • L DUT ) + exp {-j ■ J3 (v) • L) | 2 

= 2 [l + cos {P DUT (v) L DUT - p (v)L)} (3.8) 

10 

Whereas p is the phase constant of the second path, whereas Pdut is the 
phase constant of the first path through the grating 12. In this formula the fact 
is neglected that the grating 12 can have a reflection coefficient smaller than 
one. However, this has no importance for the changing part of the detected 
15 signal 36. For this changing part it can be written with (|>=p-L: 

U~ det (v)=2-cos(^ t/r (v)-^(v)) (3.9) 
Hilbert-transforming this signal gives the analytic signal: 

U~ et (v)=2-exp(^ t/r (v)-^(v)) (3.10) 

20 

and the phase information: 

& et (vO=<W(v)-<Hv) (3.11) 

To determine the time dependency of the signal when performing a frequency 
25 sweep of the laser 14 and to find an evaluation rule for the group delay T giDUT , 
the functions <)> and <)>dut are now examined. 
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For the phase of the path of the beam 28 without the grating 12 which path is 
assumed to be dispersion free it is possible to perform a Taylor approximation 
of the first order: 

</){y) =0(v o )+^-\v = v o -(v-v 0 ) 
dv 

=<t(v 0 ) + 2xjt\a> = co Q .(v-v 0 ) (312) 
=<f>(v Q ) +t(v 0 )-2x(v-v 0 ) 

=,(v 0 )3^>.2*<v-v 0 ) 
c 

where L is the length of the path without the grating 12. 

The same goes for the path 32 of the beam 26 with the grating 12: 



<W(v) =^(v 0 )+^|v=v 0 -(v-v 0 ) 



i v 



(3.13) 



2 Jv/ 

10 



where vj/ is denoting the residual error. Performing the same steps as in formula 
3.12 one reaches the following formula: 



4>DUT ( V ) = 4 DUT (Yo ) + 7 DUT Oo ) ' ^ (V - V Q ) 

+ 2^7 (W (V)) H V ° • 2;r ( v/ " v o) 2 +^ (3 ' U) 



dr DUT _dt DUT dX_ dv DUT 1 _ afr /I 2 
</v <M. 'dv dk A_(c\ dX c 

dx{x) 

20 By putting this formula into 3.1 1 this gives for the phase of the analytic signal: 
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+ In [t dut (v 0 ) - t (v 0 )] • (v - v 0 ) 
rfr. „ A 2 . ^ 2 (3-16) 
dA c 



-A^ 0 (3.17) 

+ 27rAr 0 -(v-v 0 ) (3.18) 



~\K—*<y-vtf (3-19) 



n 



+ ¥ (3.20) 

To determine the group delay of the grating 12 it is sufficient to differentiate the 
face <))ciet with respect to co. This gives: 



ftj ^del _ 1 ^Vdct _ ; r Y PUT «Y i (3 21) 

£ </<y 2;r </v 2;r r/i' </v 

O ™' 
RJ 20 

This is the information searched for since only the behavior of the group delay 
with respect to the frequency is interesting, rather than an absolute value. 

However, it is easy to extract from formulas 3.17 - 3.20 information about a 
deviation of this linear rise of the group delay with respect to the wavelength for 
25 gratings with linear chirp. This goes by approximating the measured phase 
signal with a polynomial of second order and subtracting this polynomial of the 
phase signal. After this only the part 3.20 remains which represents the 
deviation of the group delay with respect a linear chirp grating. 
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With the help of the diagram of Fig. 3 the main steps of an embodiment of the 
inventive method are described, whereas some more details of these steps 
follow later on: 

In Fig. 3 ovals describe signals and boxes describe measures of the evaluation 
unit 44. First of all the second or reference signal 56 and the first signal 36 are 
towed to the zero line in step 100 (see Fig. 4). The resulting signals 102 and 
104 are now oscillating about the zero line. In step 106 any non-linearities in 
the sweep of the laser 14 are corrected with the help of the second signal 56 
(see Fig. 5). Resulting is a signal 108 corrected for any non-linearities in the 
laser sweep. With the help of the aforementioned formulas are evaluated in 
step 1 10 the group delay 112 and the linear deviation of linear chirp 114. 

With the help of a filter step 1 1 6 it is possible to get filtered signals 1 1 8 and 1 20 
of the group delay and of the deviation from the linear chirp of the grating 12. 

In Fig. 4 the principle of the function 100 is graphically illustrated. A is the 
amplitude of the signals 36 and 56. 152 denotes the reflection spectrum of the 
grating 12. 124 denotes the ideal center line of the spectrum 122. 

The function 100 is illustrated in Fig. 8 and comprises the following steps: 

- making the signal 36, 56 oscillating about a zero line by: 

- determining the points of mean value of the signal 36, 56 by: 

- determining the maximum and the minimum of the signal 36, 56 in a 
predetermined first range of time smaller than the total range of time, 

- determining a mean value between the maximum and the minimum, 

- determining the maximum and the minimum of the signal 36, 56 in a 
predetermined next range of time adjacent the already examined range of time, 
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- determining a mean value between the maximum and the minimum, 

- repeating the last two steps until the complete time interval is covered, 

- choosing the predetermined range of time by: 

- determining the average period of the oscillations of the signal 36, 56, 

5 - choosing the size of the range so that more than two average periods fit in the 
chosen range of time. 

Figure 5 illustrates a part of the function of step 106 to correct unwitting non- 
linearities in the sweep of the laser 14. The step 106 comprises the following 
steps (in the following the symbol "f is used instead of " v "): 

10 - deriving the non-linearity information by: 

- transforming the second signal 56 to get a Fourier transformed second signal, 

- eliminating the negative parts of the Fourier transformed second signal to get 
a non-negative Fourier transformed second signal, 

- retransforming the non-negative Fourier transformed second signal to get an 
15 analytic signal of the second signal 56, 

- determining the phase of the analytic signal to get as a second phase signal 
the phase as a function of time of the second signal 56, 

- using the second phase signal for determining as the non-linearity information 
the frequency as a function of time of the second signal 56, 

20 - deriving a first phase signal by: 

- transforming the first signal 36 to get a Fourier transformed first signal, 



-15 



Agilent Technologies, Inc. 
Our Ref.: 20010487 

- eliminating the negative parts of the Fourier transformed first signal to get a 
non-negative Fourier transformed first signal, 

- retransforming the non-negative Fourier transformed first signal to get an 
analytic signal of the first signal 36, 

- determining the phase of the analytic signal to get as a first phase signal as 
a function of time of the first signal, 

- correcting the effects on the first signal 36 caused by the non-linearity by: 

- using the non-linearity information to interpolate the first phase signal of the 
first signal 36 on a linear scale of frequencies to get a corrected first phase 
signal, 

For determining the non-linearity information the following steps are performed: 

- determining the frequency f(n) of the second signal 56 as a function of n 
discrete points of time, n= 1, N, on the basis of the second phase signal by: 

- determining the second phase signal <J>(n) at the n points of time, 

- determining the maximum <j) max of the second phase signal, 

- using a predetermined maximum frequency f max of the frequency range, a 
predetermined average tuning velocity during tuning the frequency and the 
maximum (j> max of the second phase signal to determine for each of the n points 
of time the frequency f(n) according to the formula: f(n)=[(f max -f min )/(() max ]())(n), 

For the interpolating step, the following sub-steps are performed: 

- getting the linear scale f| in (n) of frequencies by: 

- using the predetermined maximum frequency f max of the frequency range and 
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the predetermined minimum frequency f min of the frequency range to determine 
the linear scale f| in (n) of frequencies according to the formula: fiin(n)=[(f max - 

fmin)/(N-1)]n, 

- preferably sorting the absolute values of f(n) monotonically, 

5 - using f(n) for interpolating the first phase signal of the first signal 36 on the 
linear scale of frequencies f| in (n). The interpolation can be done linearly or with 
the help of polynoms or splines. 

The resulting signal 108 which is corrected for any non linearities in the sweep 
of the laser 14 is used in step 1 10 for evaluate the group delay 112 and the 
M- 10 deviation 1 14 from the linear chirp of the grating 12. This evaluation is done 
with the help of the above-mentioned mathematical formulas and can 
comprises deriving a group delay 1 12 of the optical device under test 12 by 
differentiating the compensated first phase signal with respect to the frequency. 

The method can further comprise the following steps: 

[/ 15 - ignoring at the begin of the tuning a predetermined amount of values of the 
m corrected first phase signal to eliminate teething troubles out of the corrected 
% first signal 108, 

- approximating the group delay 1 1 2 with polynoms of at least second order to 
get an approximated group delay, 

20 - subtracting the approximated group delay from the group delay 1 12 to get a 
non-linear part of the group delay 1 1 2. 

- using the non-linear parts of the group delay 112 to determine the mean 
signal power of a deviation from a linear of the group delay 1 1 2 of the device 
under test 12. 
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- using the square coefficient of the polynomial to determine the mean gradient 
of the group delay 1 1 2. 

Having evaluated the group delay 112 and the deviation 1 14 it is possible but 
not necessary to make a filter step 1 16 to filter the signals 1 12, 1 14. Actually, 
this filter step is preferably performed before the calculating of the group delay 
1 12, as follows: 

- adapting the filtering to the shape of the corrected first signal 108 by: 

- a: making an interferometric signal out of the corrected first phase signal, 

- b: Fourier transforming the interferometric signal to get a spectral signal, 

- c: determining a fraction, preferred the half, of the maximum of the spectral 
signal, 

- d: determining the abscissas of the intersections of the ordinate of the fraction 
with the curve of the spectral signal, 

- e: determining the mean frequency f mean as the average of the abscissas, 

- f: band-pass filtering the spectral signal with a band-pass filter having its 
center at the mean frequency and having a width greater than the width of the 
frequency range, 

- determining the width of the band-pass filter by: 

- a: predetermining or estimating the maximal range GD ran ge of the group delay, 

- b: determining the mean value GD mean of the group delay according to the 
steps c-e of claim 26, 

- c: calculating the filter width according to the formula: filter width= 
- 18 
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f mean(G D ran ge/ G D me an )i 

- filtering with a filter having the shape of a Hanning window. 

Resulting is a filtered group delay 118 and a filtered deviation 120. 

Figure 6 shows the principle of determination of the center frequency f mean of 
the interferometric spectrum. 

Figure 7 shows a graph showing different filter shapes for use with the inventive 
method, a flat top filter shape 130 and a Hanning form filter shape 140. Figure 
9 shows the group delay 1 1 2 over wave length of a chirped Bragg grating 1 2, 
which is a FBGDC-HP1A of the company JDS Fitel. It is clearly seen that the 
graph of the group delay 1 12 is very smooth with respect to graphs determined 
with methods of the prior art. 
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